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Abstract 

The spectral statistics of low-lying states of fp shell nuclei are studied by per- 
forming large shell-model calculations with a realistic nuclear interaction. For Ca 
isotopes, we find deviations from the predictions of the random-matrix theory which 
suggest that some spherical nuclei are not as chaotic in nature as the conventional 
view assumes. 



1. Introduction 

In recent years many authors have shown great interest in the fluctuation proper- 
ties of energy levels 1 ) . It is well known that the fluctuation properties of quantum sys- 
tems with underlying classical chaotic behaviour and time-reversal symmetry agree 

1 This work has been partially supported by the Ministero dell'Universita e della Ricerca Scien- 
tifica e Tecnologica (MURST) and by the INFN-CICYT Agreement. 
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with the predictions of the Gaussian Orthogonal Ensemble (GOE) of the random- 
matrix theory, whereas quantum analogs of classically integrable systems display the 
characteristics of the Poisson statistic 2 ) . 

In atomic nuclei, the fluctuation properties of energy levels are best studied in 
the domain of neutron and proton resonances near the nucleon emission threshold, 
where a large number of levels of the same spin and parity in the same nucleus are 
present, and an excellent agreement with GOE predictions has been found 3 ). In the 
ground-state region, however, the samples of consecutive experimental levels of the 
same spin and parity in any one nucleus are quite small. Therefore it is more difficult 
to calculate reliable mean values and fluctuations of statistics such as energy level 
spacings. In order to circumvent this difficulty, in recent years 4,5 ' 6 ) statistical analyses 
of experimental lower excitation energies have combined data from a large range of 
excitation energies and angular momenta of a nucleus or a set of nuclei. Such analyses 
have provided evidence suggesting that spherical nuclei show level spectra close to 
GOE predictions and deformed nuclei show strong deviations from GOE behaviour. 
In a recent analysis 7 ) of level spacings close to the yrast line of deformed nuclei with 
Z = 62-75 and A = 155-185, the average level spacing for states with the same 
spin and parity was calculated for the total ensemble instead of for individual nuclei. 
The level spacing fluctuations obtained are quite close to the Poisson distribution, 
showing evidence of regular motion. 

We can conclude that experimental nuclear spectra suggest a dominance of chaotic 
motion in most cases and that regular motion seems to be prevalent in a small region 
of excitation energy just above the yrast line of deformed nuclei. But the borderlines 
in mass number, excitation energy, etc., between order and chaos in nuclei are not 
sharply distinguished by the experimental data available. 

2. Analysis of the fp shell 

Theoretical models should also help to establish the domains of chaos in nuclei. 
First of all, the success of the nuclear shell model and the collective models of nuclear 
structure at low excitation energies provide strong evidence of regular motion in 
nuclei. Theoretical studies of deformed nuclei have shown that a two-body interaction 
of reasonable strength causes chaotic spectra a few hundred keV above the yrast line 
due to rotational band mixing 8 ). The transition from order to chaos has also been 
studied in the framework of the Interacting Boson Model and was found to depend 
on the values of the model parameters. Regular features of the spectra were found 
to be associated with values of the parameters for which the hamiltonian displays 
a symmetry 9 ). However, statistical analyses of shell model energy spectra and wave 
functions have almost always shown that chaos in nuclei is quite widespread 10 ). 

It is somewhat surprising that shell model calculations suggest such a prevalence 
of chaos because the model itself is based on a mean field idea associated with regular 
motion. The residual interaction can destroy the regularity of single-particle orbits 
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and cause the chaotic features, but it is still remarkable that chaos can be so dominant 
even for the lowest energy levels of light nuclei 11 ^ . Recent shell-model calculations by 
Bae et al. 12 ^ have shown that for nuclei of mass A = 212 shell model spectra display 
features of regular motion in some cases. The main reason for this seems to be the 
relatively small values of the residual interaction matrix elements as compared to the 
average spacing between neighbouring single-particle levels in heavy nuclei. 

Table 1: m-scheme and maximal (J,T) dimensions of the analyzed configurations 
(n,p) of active nucleons. 





my 




46 Sc 


m Ca 


A8 Ca 




(n,p) 


(3,3) 


(4,2) 


(5,1) 


(6,0) 


(8,0) 


(10,0) 


m-scheme dimension 


121 440 


86 810 


30 042 


3 952 


12 002 


17 276 




(4,0) 


(4,1) 


(4,2) 


(4,3) 


(4,4) 


(4,5) 


(J, T) dimension 


4 750 


8 026 


3 783 


615 


1 755 


2 468 



In this work, we undertake the statistical analysis of the shell-model energy levels 
in the A = 46-50 region. Exact calculations are performed in the {pi/2,P3/2,f 5/2^7/2) 
shell-model space, assuming a m Ca inert core. The diagonalizations are performed 
in the m-scheme using a fast implementation of the Lanczos algorithm through the 
code ANTOINE 14 ). The interaction we use is a minimally modified Kuo-Brown force, 
explained in 15 ) . For a fixed number of valence protons and neutrons we calculate the 
energy spectrum for projected total angular momentum J and total isospin T. 

Table 2: Number of energy levels up to 4 and 5 MeV above the yrast line in 46 V. 



J 


4 MeV 


5 MeV 





1 


2 


1 


7 


11 


2 


5 


8 


3 


14 


23 


4 


7 


11 


5 


13 


23 


6 


8 


13 


7 


18 


33 


8 


9 


14 


9 


14 


21 


0-9 


96 


159 



We calculate the T = T z states from J = to J = 9 for all the combinations (n, p) 
of 6 active nucleons, i.e. (3,3), (4,2), (5,1) and (6,0). For the sake of completeness 
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we also calculated (8,0) and (10,0). Table 1 shows the m-scheme and the maximal 
(J,T) dimensions of the configurations (n,p) of active nucleons analyzed. 

Since we are looking for deviations from chaotic features, we are mainly interested 
in the low-lying levels, up to a few MeV above the JT yrast line. Table 2 shows the 
number of levels up to 4 and 5 MeV above the yrast line in m V . As can be seen, 
in many cases the number of levels is too small to calculate fluctuations around the 
average spacing between neighbouring levels. 

Hence in those cases we use a larger set of levels for the calculation of the mean 
level spacing as a function of energy. For each JT set of levels the spectrum is 
mapped into unfolded levels with quasi-uniform level density by using the constant 
temperature formula 5 - 1 . Under the 8 MeV range of energies the mean level density 
can be assumed to be of the form 

p(E) = ^exp[(E-E )/T], (1) 

where T and E are constants. For fitting purposes it is better to use not p(E) but 
its integral N(E). We take 

N(E) = [ E p{E')dE' + N = exp [(E - E )/T] - exp [-E /T] + N . (2) 
J o 

The constant N represents the number of levels with energies less than zero. Fol- 
lowing Shriner et al. 5 \ we consider Eq. 2 as an empirical function to fit the data and 
let N take non-zero values. The parameters T, E and N that best fit N(E) are 
obtained by minimizing the function: 

G(T, E , N ) = f EmaX [N(E) - N(E)] 2 dE, (3) 

where N(E) is the number of levels with energies less than or equal to E. The 
energies E m i n and E max are taken as the first and last energies of the level sequence. 
As an example, Fig. 1 illustrates the fit to the integrated level density N(E) for the 

jn T = g +1 levelg of 46 TL 

Once the best fit F(E) to N(E) is obtained, the unfolded energy levels are given 

by 

E- = E ■ + (E -E ■) (4) 

^mm ~ 7-i \ 771/ T7i \\ J - J max ^min)- V / 

" \ I- 1 max ) " V f-'min ) 

These transformed energies should now display on average a constant level density 

Pc- 

The spectral statistic P(s) is used to study the local fluctuations of the energy 
levels 16 ' 17 - 1 . P(s) is the distribution of nearest-neighbour spacings Sj = {Ei + \ — Ei)p c 
of the unfolded levels E{. It is obtained by accumulating the number of spacings that 
lie within the bin (s, s + As) and then normalizing P(s) to unity. 

For quantum systems whose classical analogs are integrable, P(s) is expected 
to follow the Poisson limit, i.e. P(s) = exp (— s). On the other hand, quantal 
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Figure 1: The best fit to the integrated level density for 46 Ti with J = 6 + and isospin 
T - T z = 0. 



analogs of chaotic systems exhibit the spectral properties of GOE with P(s) = 
(7r/2)sexp(-fs 2 ) 1 > 2 ). 

The P(s) distribution is compared to the Brody distribution 18 ) 

P(s,u) = a(uj + l)s w exp(-as w+1 ), (5) 

with 

a = (r^ir 1 . (6) 
■+- i 

This distribution interpolates between the Poisson distribution (uj = 0) of integrable 
systems and the GOE distribution (u — 1) of chaotic ones. 

In order to obtain a meaningful statistic, P(s) is calculated using the unfolded 
level spacings of the whole set of J = 0-9 levels for fixed T up to a given energy 
limit above the yrast line. Thus the number of spacings included is reasonably large. 
For example, up to 4 MeV there are 52 spacings for T = and 86 for T — 1, and 
up to 5 MeV there are 87 and 149 spacings, respectively. Table 3 shows the Brody 
parameter ui for the J = 0-9 set of level spacings for the A = 46 nuclei up to 4 and 
5 MeV above the yrast line. 

Clearly, m V and m Ti are highy chaotic, but there is a considerable deviation from 
GOE predictions in 46 Ca, which is a single closed shell nucleus. Similar calculations 
for 48 Ca and 50 Ca yield uj values similar to 46 Ca, as shown in Table 3. Thus, for 
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Figure 2: P(s) for low-lying levels of the fp shell with < J < 9 for (n,p) configura- 
tions: (3, 3) and (4, 2). The dotted, dashed and solid curves stand for GOE, Poisson, 
and Brody distributions, respectively. 



the Ca isotopes we find the same kind of phenomenon obtained by Bae et al. 12 ) in 
the heavy single closed nuclei 212 Rn and 212 Pb, namely that low-lying states deviate 
strongly from chaoticity toward regularity. 

Table 3: Brody parameter uj for the fp shell for different configurations (n,p) of 
active nucleons with < J < 9 and E < 4 MeV (up) and E < 5 MeV (down). 



my 


46 T , 


46 Sc 


46 Ca 


48 Ca 


50 Ca 


0.97 
0.97 


1.09 
0.97 


0.65 
0.93 


0.57 
0.57 


0.64 
0.50 


0.56 
0.65 



To obtain a better estimate of the Brody parameter, we can separately combine 
spacings of different nuclei. In fig. 2 we plot -P(s) for the 236 spacings of 46 V + A6 Ti, 
and in Fig. 3 for the 266 spacings of 46 Ca+ A8 Ca+ 50 Ca up to 5 MeV above the 
yrast lines. The number of level spacings is now sufficiently large to yield meaningful 
statistics and we see that Ca isotopes are not very chaotic at low energy, in contrast 
to other nuclei in the same region. 

How may these results be explained? We observe that the two-body matrix el- 
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Figure 3: P{s) for low-lying levels of the fp shell with < J < 9 for (n,p) configura- 
tions: (6,0), (8,0) and (10,0). The dotted, dashed and solid curves stand for GOE, 
Poisson, and Brody distributions, respectively. 



ements of the proton-neutron interaction are, on average, larger than those of the 
proton-proton and neutron-neutron interactions. Consequently the single-particle 
motion in nuclei with both protons and neutrons in the valence orbits suffers more 
disturbance and is thus more chaotic. The shift towards the Poisson distribution 
obtained for the Ca isotopes could also be due to some underlying symmetry domi- 
nating the dynamics of those nuclei. However, we do not find significant differences 
between J = and high J values when they are analyzed separately. Thus the 
seniority scheme does not seem to have a significant influence on the shift towards 
regularity observed in Ca isotopes. 

It should be noted that the analysis 19 -* of experimental energy levels below 4.3 MeV 
excitation energy in the semi-magic nucleus U6 Sn yields a near-neighbour spacing 
distribution which is intermediate between GOE and Poisson, with uj = 0.51 ± 0.19. 
This result is consistent with the theoretical findings of Bae et al. 12 ) for 212 Rn and 
212 Pb, and our present results for Ca isotopes. 

3. Conclusions 

Why do all shell model calculations give chaotic features for sd shell nuclei, with- 
out any significant deviations towards regularity? First, it should be noted that most 
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of these calculations include a large number of states, up to excitation energies far 
above the nucleon emission threshold. This is, for example, the case of the 22 cal- 
culations of Bae et al. 12 **, which include the full set of levels for several J values and 
obtain uj = 0.96. Second, we notice that Ormand and Broglia 11 - 1 obtained a GOE-like 
distribution for the first two spacings of each spectrum for a set of sd shell nuclei. 
However, these nuclei have valence protons and neutrons and are thus similar to the 
case of 46 V and 46 Ti, for which we also find chaotic behaviour. We conclude that no 
regular features have been found in the sd shell region. This is because single-closed 
nuclei have a small configuration space and too few low-lying levels for statistical 
analysis, and also because the disturbance of single-particle motion by the two-body 
interaction is greatest in light nuclei. 
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